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Abstract 

Starting from the static Fukuyama-Lee-Rice equation for a three- 
dimensional incommensurate charge density wave (CDW) in quasi one- 
dimensional conductors a solvable model for local phase pinning by im- 
purities is defined and studied. We find that average CDW energy and 
average pinning force show critical behaviour with respect to the pinning 
parameter h. Specifically the pinning force exhibits a threshold at h = 1 
with exponent (3 = 2. Our model examplifies a general concept of local im- 
purity pinning in which the force exerted by the impurity on the periodic 
CDW structure becomes multivalued and metastable states appear beyond 
a threshold. It is found that local impurity pinning becomes less effective 
at low temperatures and may eventually cease completely. These results 
are independent of spatial dimensionality as expected for local impurity 
pinning. Comparison with Larkin's model is also made. 

PACS: 71.45.Lr, 71.55.Jv, 72.15.Nj 

1 Introduction 

In charge density waves (CDW) which appear below the Peierls transition temper- 
ature in quasi one-dimensional metals pinning of the order parameter phase 
at point defects is an important effect. Phase pinning results from the electro- 
static coupling between the spatially periodic charge modulation in the CDW and 
the electric potential of an impurity. The phase ip determines the position of the 
CDW with respect to the host lattice. The energy and force contributions of an 
individual impurity become periodic functions of if. Distortions of the phase near 
the impurity produce a positive elastic energy. By properly adjusting the phase 
pattern near the impurity a net pinning force can be possible. Phase pinning in 
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CDW is not fully understood at present. Well known is the weak or collective 
pinning limit. In this case each impurity only slightly distorts the local phase, 
however, many impurities act coherently in a macroscopic metastable Lee-Rice 
domain p|, Q and produce a small pinning force. This mechanism has its ana- 
logue in pinning of flux lines in type II superconductors |[. The opposite 
limit of strong pinning is obtained when the pinning strength of a single impurity 
is so large that the elastic energy of local phase deformations is neglegible 
Intermediate concepts have also been explored @, 0, . 

In this paper we point out a number of special features for local phase pinning 
by impurities. Local phase pinning describes pinning effects linear in the con- 
centration of impurities. We study a three-dimensional incommensurate CDW 
within a solvable model. We find that the averaged CDW energy shows a weak 
singularity as function of the pinning parameter h aX h = 1 characterized by an 
exponent h which is calculated. The formal limit of strong pinning is reached 
only for very large h. The pinning force has a threshold sX h = 1, i.e., it vanishes 
for h <1 and then starts with a "critical" exponent similar to a phase transition. 
Our model realizes a general concept of local impurity pinning when the force 
exerted by the impurity on the periodic structure (CDW in our case) becomes 
multivalued and metastable states appear |p. In a single chain model a similar 
problem has been considered recently in [P, |l3l neglecting the screening. 



Larkin's model has been extended in to include dynamics. It has 
meanwhile evolved into a complete concept of CDW and SDW dynamics |jll], ^ . 
Here, we restrict ourselves to the static properties of a complementary model. 

The model is introduced in Sec. 2 . It is exactly solved for average energy 
(Sec. 3) and average pinning force (Sec. 4) in the screened limit when an abun- 
dance of quasi-particles eliminates Coulomb forces. The latter act between the 
local charges that are produced by phase deformations in CDW. In Sec. 5 we 
discuss descreening in semiconducting CDW. We show that in lowering the tem- 
perature local impurity pinning becomes less and less effective and may even 
cease completely. Finally Sec. 6 gives a detailed account of Larkin's model 
and its relation to our results. 



2 Pinning Model 

Following we write the phase dependent part of the static energy density of 
incommensurate CDW plus one impurity at position Xj as 



The stiffness constant in the fully screened limit and for a rectangular lattice is 
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where at^ are the interchain distances of the chains which run along the z- 
direction. The Fermi velocities for the anisotropic metal are denoted vp, Va 
and Vt2, respectively. In terms of the Fourier component V(Q) of the short range 
impurity potential (Q denotes the CDW nesting vector) and in terms of the half 
gap A, the electron density no, Fermi energy ep, and electron-phonon coupling 
constant \q the pinning amplitude Vq is given by Vq = nQV{Q)A/{4:XQep) [Q. 
Vq has dimension of energy. We neglect dislocation lines in the CDW lattice and 
consider the phase to be unique over the crystal. 

As in we perform an isotropy scaling in the static CDW energy E = J (PxTi 
according to z' = z, x' = xvp/va^ and y' = yvp/vt2 and use the old names. Then 
H becomes 



n = ^K{V^y + Vo5{^ - X,) cos(Qx + (3) 

with K = hvtiVt2/{27t'atiat2Vp), while Vq remains unchanged. The formal solution 
of the Poisson equation following from (^) gives divergent energy in three dimen- 
sions. We introduce a cutoff ^ which is of the order of the amplitude coherence 
length hvp/A into our model by setting y9(|x — Xj| < C,/2) = ipi. Thus the phase 
in our quasi-isotropic setting is constant inside a sphere of diameter ^ centered 
at the impurity. ^ is a constant parameter of the model. This case corresponds 
to the following modified inhomogenity in the Poisson equation associated to (^: 

-VoS{\x-^i\-C/2) . (4 

For convenience we fix the phase at infinity by requiring (poo = 0. A phase 
V^oo 7^ can be trivially transformed away from all our equations, especially from 
the energies to be calculated. The latter, therefore, do not directly refiect pinning. 
Phase pinning, however, shows up in the pinning force f^, [T3[| . 
The solution of the modified Poisson equation for the phase is 

V.(|x-x,|>e/2) = |^^. (5) 

In (|^) Qxj acts as a random phase uniformely distributed in (— vr, vr) and we will 
denote it as P. By averaging the total energy over P we obtain the CDW energy 
per pin of a random ensemble of local pins. 
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The total CDW energy associated with the solution fi^x) becomes E = T^K^ipf + 
Vq cos(r + (fi). With the pinning parameter 



the energy can be expressed as 



^ = ^° + cos(r + y,,)^ . (7) 
Finally the phase (fi follows from the requirement dEjdi^i = or 

Lfi = /isin(r + (fi). (8) 

This is an implicit equation for ipi in terms of F. Equation (H) also appears in 
the mean field model of Fisher [T3| for the multivalued static solutions a.t h > 1 



in the presence of an external field. 

Knowing (pi{T) the reduced averaged CDW energy e{h) for example can be ex- 
pressed as 



e{h)=<E> /Vo = — r ciF{-sin2(F + </.,(F)) + cos(F + (^,(F))}. (9) 

It is noted that (§) can also be obtained without considering ipi as a variational 
parameter by starting from the inhomogeneous phase equation. 
We summarize the principal features and assumptions of our model: 

1. Three-dimensional Fukuyama-Lee-Rice model with one randomly placed 
impurity 

2. Averaging over the impurity position amounts to an ensemble average over 
many local, i.e., independent pins. Weak pinning is excluded. 

3. Neglect of quantum tunneling and thermal fluctuations. For a discussion of 
the corresponding conditions see [^]. It is also noted that Vq ^ ksT holds 
for strong pins. 

The following two sections derive exact results from this model. 
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3 Average CDW Energy 



We begin by investigating the average CDW energy in the equihbrium state. This 
also serves as an introduction to our method of solution. 

Our model is nontrivial because of equation (H) and its variety of solutions. 
For h < 1 there exists one and only one solution ipi for every F. Writing (§) as 



w 



r + hsinw] w = (pi + r, 



(10) 



the solution w{r) maps the interval — tt < F < tt one to one onto — tt < u; < tt. 
It is then easy to calculate e{h < 1). Transforming the integration over F in (|^) 
to an integration over w using dT / dw = 1 — h cos w gives immediately 



1 

e{h < 1) = — / dw{l — hcosw) 

2ll J -TV 



h 2 

— sm w + cos w 



h 
I' 



(11) 



This argument fails for h > 1 when more than one solution of (|TDp exists. For 
/i > 1 at least three solutions solutions exist. For h = 3 the situation is depicted 
in Fig. 1. In general these solutions give different energies. The physical relevant 
solution is the one with the lowest energy. Choosing this solution for each F 
defines the integration path. 

Formally the average energy above the threshold can be expressed as 



e{h) 



27V 



E 



h 3h h"^ ' 

sm w w sm 2w sm w 

4 8 6 



„(■') 



(12) 



{wl'^\w^'') are the endpoints of appropriate integration intervals along the w-axis 



giving minimum energy. The endpoints of the path are w, 



(1) 



-TT, W. 



(1) = - 



Wa 



and w. 



(2) 



Wa, w^^ 



(2) 



71 where Wa = sm{wA), i.e., there is a jump from —wa to 



Wa at F = 0. 

The value /i = 1 is a singular point in the following sense: The dependence of 
e{h) + h/A on h changes taking on locally the form of a power law with exponent 
b. Eq. (|12D can be used to calculate the exponent b defined by 



e{h) + ^ ~ e^ e = h-l<^l. 



(13) 



The value of wa near h = 1 can be found perturbatively. To lowest order with 
respect to e it is wa = y/Oe. Eq. (|12D then gives 



e{h) + - 



357r 



(14) 
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i.e., the exponent is 6 = 7/2. To complete the analytical investigations, the 
behaviour for /i ^ 1 is studied. In this limit the value of wa_ is about ±7r(l — 1/h) 
which gives 



e(h) = -l+\^—-lj/h, /i>l. (15) 

The usual strong pinning limit < E >sp= —Vq is thus asymptotically approached 
albeit at a rather slow rate. It is also noted that for | arccos(l//i) — \/h'^ — 1\ > vr 
additional solutions exist which lead to more and more metastable states. 



4 Average Pinning Force 

Physically more important than the average energy e{h) is the average pinning 
force the CDW experiences when it is moved adiabatically. This force is 



The pinning force F{h) vanishes for h < 1 since dE/dV then is a single valued 
and periodic function over the interval — vr < F < tt. Physically it means that 
the pinning forces from different local impurities cancel 0. 

Above threshold transitions between branches of different energy occur. Using 
the energy profile 



e(w^,r) = — 



[w 



2h 



+ cos w 



(17) 



and 



de{w(T),T) w-T . 
^1^ = -^ = """^' 

the reduced force f{h) can be treated in anology to the energy. The result is 



1 

2^ 



'cos2wb — cos2wc) — ( cos wb — cos Wc 



(19) 



Here wb is the point of the energy minimum e{wB, Pc) into which a a transition 
occurs from the critical metastable state characterized by a horizontal inflection 
point in the energy profile at (ti^cPc)- The r.h.s. of (|19D is just the difference in 
energy of these two states devided by 27i. 
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It is possible to give a more geometrical description of this process and the 
resulting force: Above threshold the force from shifts of the CDW by less than 
X = Tc/Q where Tc is given by Tc = ^h? — \ — Wc with Wc = arccos(l//;.)) is 
elastic. Larger shifts populate metastable states. 

For Fc < vr the averaging path is from — vr to —Wc and then from wb to tt as 
shown in Fig. 1. It runs over the bifurcation point (F^ —Wc) and jumps to the 
upper branch starting at wb = Fc + hsm{wB)- This is clearly seen in the energy 
profile e{w,T). 

Near h = 1 one can use the perturbative results Wc = wa = '^Wc to find 



F{h) = ^QV^e\ e«l. (20) 
47r 

Thus pinning sets in at the threshold h = 1 with exponent f3 = 2. 

For larger h the situation is more complicated since more metastable states 
exist. The averaging path now extends over n periods in order to sample all 
metastable states. The jump is from (F^, —Wc) to the point (F^ ws) which gives 
the lowest energy, wb is that value on the curve F + hsm{w) that hes nearest 
to (2n — 1)tt where E{w) has its absolute minimum. The energy difference is 
E{-Wc) - E{wb) with E{±Wc) = Voih/2 + l/(2/i)). Whenever F^ crosses a 



value 27m {n = 1, 2, 3, ...), i.e., when h = hn where /i„ is found from ^K^ — 1 — 
arccos(l//i„) = 2n7r (approximate solution hn = 2n7r + 7r/2) the winding number 
changes from n to n + 1 and f{h) is reduced correspondingly. This leads to a 
saw-toothed structure of f{hn-i < h < = {h/2 + l/{2h) — E{wb)) I {2Tm) 
[Hq = 1). Thus f{h) oscillates around the asymptotic value 1/2 with approximate 
period 27r and decreasing peak to peak amplitude 5/(/i„) = f{hn)/{n + 1). Fig. 
2 displays the reduced pinning force f{h). It is noted that the unsealed force 
(|T6|) contains the pinning amplitude Vq and does not saturate for — oo as in 
Larkin's model 0. 

The threshold behaviour of the pinning force is also found in the one-dimensional 
pinning model in ^ and is a general feature of local impurity phase pinning 

i- 



4.1 Pinning Parameters 

The relation of the pinning parameter h to the usual measure of pinning strength 
will now be discussed. According to (0) and taking the isotropy scaling into 
account, h is given by: 



h = Vo 



atiat2VF 

VtlVt2H ' 



(21) 



The standard measure |T^ for impurity energy to elastic energy in an anisotropic 
three-dimensional CDW continuum is 
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ei = 2W^anancT (22) 
\Vt\Vt2) nvp 

where C3 is the impurity concentration. The effective distance between pins is 
i = {c3VtiVt2/vp)''^^'^. Thus ei is related to h by 

ei = 2Txi hji. (23) 

The apphcabihty of the theory requires ^ to be smaller than I. Using the values 
fi? = 5 ■ 10 ''cms"^ ~ 10ft, A = 1000/csT, and an impurity concentration of one 
ppm the ratio ^jl is estimated as 0.01 One can, nevertheless, conclude that the 
strong pinning limit in our model, /i ^ 1, also requires large e^. Weak pinning 
occurs for < 1. The region < /i < 1 where local pinning is absent thus is 
concealed by weak pinning which is ubiquitous in less than four dimensions. 



5 Descreening 

So far we have considered the case of full screening: There are enough thermally 
excited quasiparticles (or normal carriers as in NbSea) to completely screen out 
the Coulomb forces between charge fluctuations associated with phase deforma- 
tions. It is known that descreening stiffens the CDW. This stiffening of the CDW 
leads to a corresponding increase of the phason velocity which has been observed 
by neutron scattering and explained in |TB[ as a descreening effect. In a sim- 



ple approximation involving only the condensate fraction < 1 one can define 
an effective stiffness constant 



K 



(24) 



T9|,pO|,P, 123,110 



to take care of descreening within the elastic CDW model 
[2^ . Following ||2^ this is a reasonable approximation for 1 — > C = Q^^f/ (Scq) 
[et- static transverse dielectric constant). In the opposite limit Coulomb inter- 
actions require a different approach. The modified stiffness constant changes the 
pinning parameter h according to 



h^hVl-N = heff, (25) 

since a scaling in chain direction becomes necessary to maintain quasi isotropy. 
Our earlier formulae hold with this replacement. When h^ff becomes less than 
unity in decreasing the temperature local impurity phase pinning stops. The 
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singularity in the local pinning force at heff = 1 is likely to be masked by thermal 
fluctuations and by weak pinning. 

The fact that local pinning centers require a minimum strength to act as 
strong pins is possibly the reason for the following observation in p5|: The CDW 
in Ti doped NbSes which is fully screened due to a partially gapped Fermi surface 
exhibits weak instead of the expected strong pinning. 



6 Relation to Larkins's model 

We want to point out the similarities and differences in our approach to Larkin's 
model iOH and fill in some additional information about it. 



Larkin uses a single chain model |26] with one impurity at x = Xi. It is defined 
by the energy functional: 



E 



E.. 



dx 



1 — cos(p{x) + 



1 / dip{x) 



dx 



+ Vo[l-cos{T + ipixi))]. (26) 



In ( ]26|) the length x is scaled to hvp/{2t±). The interchain coupling energy is 
denoted by t±_ and Es = Stx/vr is the energy of a 27i phase soliton {Es is 
called w m B). 



/^From the general theory in p4| which describes the CDW as a system of coupled 
chains Larkin's model follows under three assumptions: 

1. Phases on neighbouring chains are set to zero which turns the interchain 
coupling into the Sine-Gordon type self interaction 1 — cos<y9 in (|26|). Thus 
phases on different chains are independent. In contrast the phase varies 
slowly across the chains in our model, a case more appropriate to a screened 
situation. 

2. The low temperature or descreened limit is understood when the quasi- 
particle fraction 1 — is smaller than the Coulomb coupling constant (. 

3. The value ( of the Coulomb coupling constant mentioned in Sec. 5 is fixed 
at 1/8. In reality a smaller value holds because of vp = O(10^cms~^). 



The model ( PBD has uncharged dipole solutions for which the phase cpi = (p{xi) 
at the impurity obeys the matching condition {h = 4Vq/Es) 



2 sin — = h sin w. 
2 



(27) 



usually inequivalent-solution Lpf^ with ipY' = 2tt — Lpf . The elastic energies are 



This equation replaces our equation (|l3)- For any solution ip^f' it has another- 



(2) 
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Es[l - cos(<^l 72)] and E,[l - cos((^r72)] = E,[l + cos((^l 72)]. We consider 
here the domains < w < 27r and < (fi < 2tt. The model is thus characterized 
by one metastable state. 

The mechanism of pinning force generation for h < 2 is precisely the same 
as discussed in Sec. 4: Transition from a metastable state with relative energy 
AE which becomes a horizontal inflection point in the energy profile E[w, T) for 
r = Fc down to the ground state. The average pinning force is then: 



F = (28) 

In our model this mechanism prevails for all /i > 1 and more and more metastable 
states appear for increasing h. In Larkin's model no horizontal inflections points 
exist for h > 2. Instead two 27r-solitons are created when ipi changes by 27r in 
a corresponding change of F. This leads to F{h > 2) = QEs/tt in while the 
average value of F{h) increases linearily with h in our model as implied by the 
strong pinning concept. 

It is possible to study the special case h = 2 analytically because the exact 
solution of (^) - expressed as w = w(T) - is available: 

F 47r 

w = — I n, n = 0, 1; w = — F + 2tt. 

3 3 

From the energy profile 



e{w, F) = — = 1 - cos ^— + - sm^ -, (29) 



and the matching condition (|27|) one finds 



de{w{T),T) h . 1 . w-V 

; = — smw = — sm . 30 

ciF 4 2 2 ^ ' 

The inflection points for h = 2 are (F^ = 3tt/2,Wc = tt/2), (F^ = tt/2,Wc = 37r/2) 
and a corresponding ground state is (Fc = 37r/2,WB = ll7r/6). Thus one finds 



F{h = 2-0) = ^ [e{w,, Fe) - e{wB, F^)] = ^ cos ^ = ^QE,. (31) 

/TT 47r O OTT 

For h = 2 + the path to follow in the integral 



Jo dr Stc Jo ^ ' ^ ' 



27r Jo 
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is continuous and goes from w{0) = 2tt to w(Tc) = 37r/2 and then back to 
w{27i) = 271 resulting in a phase change of Aipi = —271. The integration gives 



= 2 + 0) = (33) 

TT 

This is the value of the force for aX\ h > 2 as pointed out in 0. 

Further analytical results are found for h < 2. From the topological condition 
dr/dw = of the critical inflection point one gets 

e{w^,T^) = 1 + ^-1^1^, (34) 

with 



A-h U-h n^-l 

Wc = arccos V ' Fc = arccos W + 2 arccos W — - — . (35) 

The perturbative expansion of the force F{h) slightly below h = 2 then leads to 

F{h) = F{h = 2-0) -^i^V3 + ^^^V2^ + 0{2-h). (36) 

It is clear that the threshold for pinning is h = 1 because there are no metastable 
states for h < 1. The intermediate regime 1 < /i < 2 is treated numerically 
and the result is shown in Fig. 3. The behaviour near /i = 1 is again F{h) ~ 
{h — ly = e"^. The critical exponent is thus {3 = 2 which seems to be universal 
for local pinning. 



7 Discussion 

From its very definition local impurity pinning is expected to be independent of 
spatial dimensionality d. This is born out by our approach. Repeating the cal- 
culations of Sec. 2 for d = 2 and d = 1 always leads to the central equation ([lOD . 
However, the pinning parameter h = hd is not any more given by (|^) which refers 
to d = 3. Using half the the mean distance id'^ C, between impurities as the dis- 
tance from the impurity where ip vanishes it is found that h2 = Voat hi{i2/^)/{hvt) 
for d = 2. For d = 1 hi = 7rVo-^i/(2ftVi?) is obtained. The relation (p3D between 
the pinning parameter hd and becomes ej = 27r/i2/ ln(^2/0 d = 2 and 
Si = 4 hi for d = 1. These relations are similar to ( pBj) and do not change the 
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conclusions significantly. The case d = 1 which does not require the short distance 
cut-off ^ is unrealistic since real CDW are never one-dimensional. 

In summary we have studied a solvable model of local impurity phase pinning 
which realizes the local pinning scenario in [5], namely singular points, thresh- 
old behaviour, and metastable states. In contrast to pinning in type II super- 
conductors static descreening and the possible deactivation of local pins at low 
temperatures are unique to semiconducting CDW (and spin density) systems. 
These results have been obtained within a phase only model plus some amend- 
ments for descreening. Especially at low temperatures more general models, e.g., 
those in ||2^ which take nonlinear screening (band bending) into account may be 
considered. 



The authors thank P.B. Littlewood who mitiatcd this study. They also thank S.N. Artemenko 
for helpful discussions. 
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Figure 1: Display of the solutions of the self consistency equation ([T0|) in the 
form w = r + hsinw for h = 3, i.e., not too far above the pinning threshold 
h = 1. The circles 1, 2, 3 indicate the solutions for — < F < Fc. (Fc, —Wc) 
is bifurcation point. The pinning force results from the vertical transition from 
—Wc to Wb on the upper branch. 
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Figure 2: Reduced pinning force / in units of QVq as function of pinning param- 
eter h. The local impurity phase pinning force / sets in at /i = 1 with exponent 
P — 2. Saturation value f{h — > oo) = 1/2, the quadratic treshold behaviour, 
and the jumps at the points /i„ ~ 27m + 7r/2, n — 1,2, ... are indicated by broken 
hnes. 




Figure 3: Reduced pinning force f{h) in units of QEg as function of pinning 
parameter h for Larkin's model. There is no pinning force below h — 1 and 
f{h) jumps from 3\/3/(87r) to the final value l/n a,t h — 2. Note the different 
normalization of f{h) in comparison to Fig. 2. 
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